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Abstract 

We present first results for the 2 —> 4 single-parton scattering gg cccc subprocess for the 
first time fully within the fcf-factorization approach. In this calculation we have used the Kimber- 
Martin-Ryskin unintegrated gluon distribution which effectively includes some class of higher- 
order gluon emissions, and an off-shell matrix element squared calculated using recently devel¬ 
oped techniques. The results are compared with our earlier result obtained within the collinear- 
factorization approach. Only slightly larger cross sections are obtained than in the case of fhe 
collinear approach. Inclusion of transverse momenta of gluons entering the hard process leads 
to a much stronger azimuthal decorrelation between cc and cc than in the collinear-factorization 
approach. A comparison to predictions of double parton scattering (DPS) resulfs and fhe LHCb 
dafa sfrongly suggests that the assumption of two fully independent DPS (gg —>• cc O —)■ cc) 
may be too approximate. 
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I. INTRODUCTION 


At high energy, gluon-gluon fusion becomes the dominant mechanism of heavy cc or 
bb pair production. The cross section for single pair production can be calculated either 
in collinear next-to-leading order approach or the /cf-factorization approach. The QQ and 
Higgs production are golden reactions for applications of the fcrfactorization approach 


ifll 


1411 . In the fcrfactorization approach the basic ingredients are so-called unintegrated 
gluon distribution functions (UGDFs) and off-shell mafrix elements. Different models 
of UGDFs have been proposed in the literature. The Kimber-Martin-Ryskin (KMR) lIlSll 
UGDF is believed to include the dominant higher-order corrections. The off-shell mafrix 
elements for gg —> QQ were calculated already long ago |lll-l3(]. The fcrfactorization for¬ 
malism was applied recently in the context of experimental data measured at the LHG 
ll6l - |20tl and a relatively good description was obtained when using the KMR UGDF. 

In the case of the Higgs boso n p roduction both 2—^1 and 2—^2 subprocess have to be 
taken into account ||14|| . In Ref. |2l|[ ] a 2 —)• 3 gg —)■ ccy subprocess was taken into account 
when calculating cross sections for pp —)■ cc'yX reaction. Recently the fcf-factorization 
approach was also applied to three-jet 11221] and Zbb 11231] production. 

A convenient formalism for the automation of the calculation of free-level scatter¬ 
ing amplitudes with off-shell gluons for arbitrary processes was recently introduced in 
Ref. |24ll . Off-shell gluons are replaced by eikonal quark-antiquark pairs, and the am¬ 
plitude can be calculated with the help of standard local Feynman rules, including the 
eikonal gluon-quark-antiquark vertex and the eikonal quark-antiquark propagator. The 
well-known successful recursive methods to calculate tree-level amplitudes can directly 
be applied, including the "on-shell" recursion, or Britto-Gachazo-Feng-Witten recursion, 
as shown in Ref. ll25l] . The heuristic introduction of the formalism in Ref. ll24|] has be giyen 


solid ground in Ref. 1261] . Most of the afford was devoted to dijet production 1271.1281] so 
far. 


The pp —?■ ccccX reaction is interesting by itself. It was shown by us recently t 


tia 


reaction is a golden reaction to study double-parton scattering (DPS) processes 129 


this 


300. 


The LHCb collaboration confirmed the theoretical predictions and obtained a large cross 
section for production of two mesons, both containing c quarks or both containing c an¬ 
tiquarks Bin . The single-parton scattering (SPS) contribution was discussed in Refs. 11321] 
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and 11331] ■ In the first case 11321] a high-energy approximation was used neglecting some 
unimportant at high energies Feynman diagrams. Last year we have calculated the 
lowest-order SPS cross section(s) including a complete set of Feynman diagrams 11331] in 
the collinear-factorization approach. The final result was only slightly different than that 
obtained in the high-energy approximation. 

In the present letter we wish to go one step further and try to calculate the SPS cross 
sections for the pp —> ccccX reaction consistently in the fcrfactorization approach. Doing 
so we may hope that a sizeable part of higher-order corrections will be included. On the 
technical side this will be a first calculation within the fcf-factorization approach based on 
a 2 —7- 4 subprocesses with two off-shell initial-state partons (gluons). The result is also 
important in the context of studying DPS as the considered SPS mechanism constitutes 
an irreducible background, and its estimation is therefore of prior importance if deeper 
conclusions concerning DPS can be drawn from measurements at the LFIC. 


II. FORMALISM 


PB 

FIG. 1: A diagrammatic representation of the considered mechanism of cccc final-state production 
via single-parton scattering within fcf-factorization approach. 

Within the fcf-factorization approach the SPS cross section for pp —)■ ccccX reaction, 
sketched in Fig. [1} can be written as 

dCpp^cccc — J dx\ dX2 H )d(Tgg^cccc ■ ( 2 . 1 ) 

In the formula above T{x,k'^,}i^) are unintegrated gluon distributions that depend on 
longitudinal momentum fraction x, transverse momentum squared of the gluons en¬ 
tering the hard process, and in general also on a (factorization) scale of the hard process 
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FIG. 2: Momenta of the off-shell gluons, represented as double lines on the left hand side, and the 
eikonal quark-antiquark pairs. The amplitude is independent of a simultaneous momentum shift 
^ibgn + <?/ ^lend + <? as long as pA'C] = 0. The same holds for the other eikonal line with ps . 



FIG. 3: Some terms in the expansion of the amplitude in terms of the eikonal propagators. The 
eikonal quarks are denoted by lines without arrows. The double lines on the left hand side rep¬ 
resent the off-shell gluons. This expansion does not represent the organization of the calculation, 
and only gives an impression of which graphs are included. 


The elementary cross section in Eq. (12.11) can be written somewhat formally as: 


da = 


d^pi d^p2 d?p2, d^p4^ 

2Ei(27r)3 2E2(27r)3 2E3{2n)^ 2E^{27i)^ 
1 


{2n)^d^{pi + p2 + P3 + PA-h- h) 


X 


flux 


\-^g*g*^cccc 


ccc (h.h)\ 


( 2 . 2 ) 


where only dependence of the matrix element on four-vectors of gluons ki and ^2 is made 
explicit. In general all four-momenta associated with partonic legs enter. The matrix 
element takes into account that both gluons entering the hard process are oft-shell with 
virtualities = —k\^ and k\ = —k^^. The matrix element squared is rather complicated 
and explicit formula will be not given here. 


As mentioned in the introduction, the scattering amp 


gluons are constructed using the formalism of Ref. 11241] . in which off-shell gluons are 


itudes with off-shell initial state 


represented by eikonal quark-antiquark pairs in order to arrive at gauge invariant ampli- 
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tudes. Figure |3] gives an idea of what kind of graphs are included. The Feynman rules 
related to the eikonal quark-antiquark-gluon vertex and eikonal propagator are 


<Wm''}i,b = 


K 


] 


Pa-K 


(2.3) 


where pA is the longitudinal momentum associated with the off-shell gluon. The external 
eikonal quark-antiquark pairs carry fundamental color indices, say i,j. It was noted in 
Ref. ll34n that the amplitude is traceless with respect to these indices, so an adjoint color 
index can be assigned to the off-shell gluon by contracting the amplitude with 
The squared amplitude summed over colors gives the same result. Denoting hy the 
amplitude with the color of one off-shell gluon highlited explicitly we have 


UMy=i: 




h] 


= E Vh ) = El-^i'' 




(2.4) 


h] 


The first term on the right hand side of Fig. |3] contains the "actual off-shell gluons" as 
virtual gluons, representing complete propagators. This term would diverge if ^ 0 
and/or —> 0, so the whole amplitude has to be multiplied with to reproduce 

correct collinear limit. 

The calculation has been performed with the help of A Very Handy LIBrary |^. In 
this Fortran library, scattering amplitudes are calculated numerically as a function of the 
external four-momenta via Dyson-Schwinger recursion 1351] . It is a recursion of off-shell 
currents, which automatically factorizes the calculation of the sum of all Feynman graphs 
such that the multiplications represented by vertices are executed only once for each ver¬ 
tex, while such vertex may occur in several graphs, for identical kinematics. This re¬ 
cursion is sketched in Fig. |4] and Fig. |5l The auxiliary eikonal quarks and anti-quarks 
are treated as external particles, so eventually an eight-point amplitudes are calculated. 
AVHLIB allows for various choices of the representation of the external helicities and col- 

nn, 

ors. These include both the color-ordered representation II361.13711 , with exact summation 
over color, and the color-dressed representation 11381.1391] . with Monte Carlo summation 
for large multiplicities. Helicity configurations can be summed exactly, or, again for large 
multiplicities, treated in a Monte Carlo approach, both discrete and with continuous ran¬ 
dom polarizations 101. The library includes a full Monte Carlo program with an adaptive 
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FIG. 4: Dyson-Schwinger recursion for off-shell currents. The thick lines represent off-shell (vir¬ 
tual) particles and the thin lines represent on-shell external particles. The sum is over all partitions 
of these external particles over the different blobs and all flavors for the virtual particles that are 
allowed according to the Feynman rules. 



FIG. 5: An explicit example of one Dyson-Schwinger recursive step for a certain off-shell current. 


phase space generator [14 ll. 14211 that deals with the integration variables related to both the 


initial-state momenta and the final-state momenta. 

The program can also conveniently generate a file of unweighted events, which ap¬ 
proach was used for the analysis presented in this paper. In the present calculation we 
use: = (Ef the factorization scale and nic = 1.5 GeV in both /cf-factorization 

and in the reference collinear-factorization calculations. Uncertainties related to the 


choice of the parameters were discussed e.g. in Ref. 11331] and will be not considered 


here. Here we wish to concentrate on the relative effect and modifications with respect 
to the results of the collinear-factorization calculations presented already in the literature 

10, S- 
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III. FIRST RESULTS 


In this section we wish to compare the new results of the fcrfactorization approach to 
those obtained by us in Ref. 11331] in the collmear-factorization approach. 

In Fig. we show standard single particle distributions in charm quark/antiquark 
transverse momentum (left panel) and its rapidity (right panel). We predict an enhance¬ 
ment of the cross section at large transverse momenta of c or c compared fo fhe collmear- 
factorization approach. The rapidity distributions in both approaches are rather similar 
(see the left panel of the figure). 




FIG. 6: Distributions in c quark (c antiquark) transverse momentum (left panel) and rapidity 
(right panel). The fcf-factorization result (solid line) is compared with the collinear-factorization 
result (dashed line). 


Distributions in rapidity of fhe cc (or cc) and cc, defined as Ycc = (i/c + J/c) /2 and Y/c = 
(i/c + i/c) /2 respectively, are shown in Fig. [71 The distributions are much narrower than 
those for single quark/anfiquark which reflecfs fhe facf fhaf fhe fwo differenf c quarks (or 
two different c antiquarks) have typically different rapidities. The discussed distributions 
in i/c arid Ycc would be identical only if i/c,l = ]/c,2 (sfrong rapidify correlations). This will 
become clearer when inspecting rapidity difference in fhe nexf plofs. 

Similar distributions but for rapidify distance between two c quarks (or two c anti¬ 
quarks) and between c and c are shown in Fig. |8l On average the distance between c 
and c is larger than that for c and c. This can be understood easily in the high-energy 
approximation discussed in Ref. ll32i] by inspecfing the contributing diagrams. Some en¬ 
hancement at small rapidity separations can be observed in the fcf-factorization approach 
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FIG. 7: Distributions in rapidity Ycc = (j/c + J/c) /2 (left panel) and Ycc = (i/c + J/j) /2 (right panel). 
The meaning of the curves is the same as in Fig. 0 

compared to the colltnear approach. 




FIG. 8: Distributions in the difference of rapidities AY^c = i/c — J/c (left panel) and AY^c =yc — yc 
(right panel). The meaning of the curves is the same as in Fig.|6l 

The distributions in rapidity distance are strongly correlated with Mcc or Mcc distri¬ 
butions shown in Fig. |9l Those distribution are, however, difficult to measure as rather 
mesons are measured and not quarks or antiquarks. 

Quite interesting are azimuthal angle correlations between c and c or c and c. The 
corresponding distributions are shown in Fig. [TOl We note much bigger decorrelation 
of two c quarks or c and c in the kf-tactorization approach compared to the collinear 
approach. This is due to explict account of gluon virtualities (transverse momenta). We 
will return to this point when discussing azimuthal correlations between mesons at the 
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FIG. 9: Invariant mass distrtributions in Mcc (left panel) and Mcc (right panel). The meaning of 
the curves is the same as in Fig. [6l 

end of this section. 




FIG. 10: Azimuthal angle correlations between two c quarks (left panel) and between c and c 
(right panel). The meaning of the curves is the same as in Fig.|6l 


Next we wish to visualize the regions of the transverse momenta of initial gluons that 
give sizeable contribution to the SPS pp —)• ccccX cross section. In Fig.[TT]we show a two- 
dimensional distribution in intial-gluon transverse momenta. The dependence on kn and 
k 2 t shown in the figure is determineated by the UGDF used in the calculation as well as by 
the dependence of the matrix element on kit arid ^ 21 - Other models of unintegrated gluon 
distributions would give different dependencies. Clearly we get large contributions from 
the regions far from the collinear case (kit = 0 and k 2 t = 0). This has of course consequences 
for other observables discussed above through the dependence of the matrix element on 
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\fs = 7 TeV 


pp^ccccX 



FIG. 11: Two-dimensional distribution in transverse momenta of initial gluons in the pp —)■ cccc 
SPS process at ^/s = 7 TeV. 


the gluon transverse momenta \Ai{kn,k 2 t)\'^ and its correlation with other kinematical 
variables. 

We will not discuss in the present letter the correlations between the gluon virtuali¬ 
ties (or their transverse momenta) and other kinematical variables related to the charm 
quarks and antiquarks. 

So far we have considered production of cccc quarks/antiquarks. As discussed in our 
previous paper 11331] such a final states may lead to the production of two D mesons, both 
containing c quarks or both containing c antiquarks which is not possible e.g. for the 
cc final-state case. As explained in Ref. |l33[l the DPS gives cross sections very similar to 
those measured by the LHCb collaboration [13ll] . How important is the SPS contribution 
discussed in this paper, calculated here in the fcf-tactorization, is shown in Fig. [121 For 


comparison we show also SPS results calculated in collinear-factorization approach liSSll . 


The two approaches give somewhat different shapes of correlation observables, inspite 
that the integrated cross sections are rather similar as discussed already at the parton 
level. Our results, so far the most advanced in the literature as for as the SPS contri¬ 
bution is considered, are not able to explain discrepancy between DPS contribution and 
the LHCb experimental data. If the discrepancies are due to simplifications in the treat- 
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FIG. 12: Distributions in invariant mass (left) and in azimuthal angle between both D^'s 

(right) within the LHCb acceptance. The DPS contribution (dashed line) is compared with the SPS 
one calculated within the Icrfactorization approach (dashed-dotted line). The SPS result from our 
previous studies 1^ . calculated in the LO collinear-factorization approach, is also shown here 
(dotted line). 


ment of DPS requires further studies including for example spin and flavour correlations. 
Some works in this direction already started 143(1 . 


IV. CONCLUSIONS 

In the present paper we have made a first calculation of the cross section for pp — 
ccccX in the fcf-tactorization approach, i.e. focussing on single parton scattering process. 
This is a first 2 —> 4 process for which krfactorization is applied. In this calculation 
we have used the Kimber-Martin-Ryskin unintegrated gluon distribution(s) which ef¬ 
fectively takes into account the dominant higher-order corrections. The off-shell matrix 
element was calculated using a new technique developed recently in Krakow. 

The results of the fcrtactorization approach were compared with the results of the 
collinear-factorization approach. In general, the kf-tactorization results are only slightly 
bigger than those for collinear approach. An exception is the transverse momentum dis¬ 
tribution above 10 GeV where a sizeable enhancement has been observed. Inclusion of 
gluon virtualities leads to a decorrelation in azimuthal angle between c and c or c and c. 

Since the cross section is in general very similar as for the collinear-factorization ap¬ 
proach we conclude that the cccc final state at the LHC energies is dominantly produced 
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by the double parton scattering as discussed in our recent papers, and the SPS contri¬ 
bution, although interesting by itself, is rather small. A comparison to predictions of 
double-parton scattering results and recent LHCb data for azimuthal angle correlations 
between and or and mesons strongly suggests that the assumption of two 
fully independent DPS {gg —)■ cc (D gg —t cc) may be too approximate or even not valid. 
Some possible reasons were discussed in Ref. 14311 . The effect found there is, however, 
too small to explain a rather large effect observed by the LHCb collaboration. This re¬ 
mains a challenge for future theoretical studies and should be confirmed by the LHCb 
collaboration at = 13,14 TeV. 
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